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SPECTRAL GAP AND EDGE EXCITATIONS 
OF d-DIMENSIONAL PVBS MODELS ON HALF-SPACES 

MICHAEL BISHOP, BRUNO NACHTERGAELE, AND AMANDA YOUNG 


Abstract. We analyze a class of quantum spin models defined on half-spaces in the d-dimensional hyper- 
cubic lattice bounded by a hyperplane with inward unit normal vector m G R'*. The family of models was 
previously introduced as the single species Product Vacua with Boundary States (PVBS) model, which is a 
spin-1/2 model with a XXZ-type nearest neighbor interactions depending on parameters Xj G (0, oo), one 
for each coordinate direction. For any given values of the parameters, we prove an upper bound for the 
spectral gap above the unique ground state of these models, which vanishes for exactly one direction of the 
normal vector m. For all other choices of m we derive a positive lower bound of the spectral gap, except for 
the case Ai = ■ • ■ = = 1, which is known to have gapless excitations in the bulk. 


1. Introduction 

One of the essential properties to understand the low-temperature behavior of a quantum lattice model is 
the presence or absence of gapless excitations above the ground state or, equivalently, whether or not there is 
a non-vanishing spectral gap above the ground state. Even when the ground state is known, answering this 
question is, in general, quite non-trivial, especially in higher dimensions. We are interested in developing 
techniques to prove lower bounds for the spectral that work in two and more dimensions by studying specific 
models. 

The issue of the existence of a spectral gap may be further complicated by the presence of gapless edge 
states, which can occur for certain geometries while the excitations in the bulk remain gapped. Edge 
states play a central role in characterizing quantum many-body states. The occurrence of low-energy, often 
gapless, states supported near the boundary of an extended many-body system have been connected with 
entanglement properties and topological order m and with phenomena such as the quantum Hall effect and 
the spin Hall effect. They may reflect the correlation structure of the state in the bulk and have a direct 
bearing on the classifcation of gapped ground state phases and the phase transitions between them [sjiniis]. 

Progress in the classification of ground state phases and in our understanding of topologically order in 
many-body models has mostly come from the study of classes of models with simplifying features. E.g, 
Kitaev’s Toric Code model [5] and the Levin-Wen models [5] are frustration-free and have the additional 
property that the Hamiltonian is a sum of commuting terms. Models with Matrix Product Ground States 
(MPS) in one dimension [3 [16] or Projected Entangled Pair States (PEPS) in two and more dimensions 
[HKniii] have also been studied with considerable success. Isolating the features of interest in models with 
ground states that are as simple as possible has proved to be a productive strategy. The present work follows 
the same philosophy. 

The Product Vacua with Boundary States (PVBS) models, introduced in one dimension in [I], were 
generalized to d dimensions in |2]. In the latter work the authors proved in a particular example that a 
model that has a non-vanishing spectral gap above the ground state in the thermodynamic limit and when 
defined on finite rectangular boxes, may have a spectral gap that tends to zero when defined on a sequence 
of diamond-shape finite volumes (with edges at 45 degree angles between the edges and the coordinate axes) 
of increasing size. It was shown that the GNS Hamiltonian for the model on such an infinite half-space is 
gapless due to edge excitations while it remains gapped in the bulk. 

In this work we extend the results of [2] by considering the spin-1/2 model (corresponding to a single 
species of particles) on half-spaces in bounded by an arbitrary hyperplane. The Hamiltonian has XXZ- 
type nearest neighbor interactions depending on parameters Xj G (0,c»), one for each coordinate direction, 
see (13) for the definition. For any given values of the parameters, we prove an upper bound for the spectral 
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gap above the unique ground state of these models, which vanishes for exactly one direction of the normal 
vector m. For all other choices of m we derive a positive lower bound of the spectral gap except for the case 
Ai = • • • = Ad = 1, which is known to have gapless excitations in the bulk. The gapless situation, while 
occurring only for one specihc orientation of bounding hyperplane, is nevertheless of particular interest. It is 
possible that an entropic selection effect occurs that makes interfaces and free surfaces in the gapless direction 
more common than one might a priori expect. Regardless, the mathematical and physical properties of the 
gapless edge spectrum deserve to be explored more fully in future work. For instance, PVBS models with 
gapless edges states in the presence of disorder would provide an interesting framework to study Many-Body 
Localization effects at interfaces. 


2. Definition of the Model and Main Results 


2.1. The single species PVBS Model. In this paper we consider the single species PVBS model on 
as introduced in [5], and adopt the notations of that paper. Let ei,..., denote the canonical basis vectors 
of and, by the natural embedding, also of For a finite lattice A C Z'^, we associate with each site 
X G A the two-dimensional Hilbert space TLx = C^, with the orthonormal basis vectors |0) and |1) describing 
a site that is either empty or occupied by a particle. The Hilbert space for the system on A is given by the 
tensor product "Ha = (Six^a'Hx- The Hamiltonian is a sum of projections hx,x+ej ,with j = 1,.. ., d, such 
that x^x + Cj G A. The nearest neighbor interactions act non-trivially only on two copies of and depend 
on parameters Xj G (0, oo), j = 1,..., d. They are dehned by 

(1) = |</.("^))(</.("^)| + |ll)(ll|, 

where = (|01) — A|10))/Vl + A^, for A G (0, oo). Here |01) is shorthand for |0) 0 |1) where the Hrst 
tensor factor is associated with the site x and the second with x + ej, etc. The Hamiltonian is defined by 


( 2 ) 


d 

^ ^ hx^x+ej 

s.t. x+gj-eA 


which is frustration-free and translation invariant. Alternative expressions for this Hamiltonian in terms of 
spin matrices and hard-core Boson creation and annihilation operators are given in [2]. 

As shown in [2], for a bounded and connected volume A, this Hamiltonian has a two-dimensional ground 
state space which is spanned by the zero particle state and a one-particle state given as follows: 


( 3 ) 


^0 = 10)"^ := Oo^gaIO) , 


1 

7W) 


xGA 


where cr^ is the first Pauli matrix acting on the site a;, and (T^|0)^ is the state with one particle at site x. 
Here and in the rest of the paper, A“ = XjA C'(A) is the normalization factor so that ||'I'j'^|| = 1. It is 

given by 


( 4 ) 


c(A) = Ea^"- 


We are interested in the excitation spectrum of this model defined on infinite half spaces bounded by a 
hyperplane containing the origin, that is subsets £> C Z"^ determined by a unit vector m G (the inward 
normal) as follows: 

(5) D := {x G : m ■ X > 0}. 


More precisely, we are interested in the spectrum of the Hamiltonian in the GNS representation of the 
zero particle ground state on D. Since this ground state is given by a tensor product vector, the GNS 
representation can be given explicitly. The GNS Hilbert space is generated by all states with only a hnite 
number of occupied sites. The dense subspace spanned by all such vectors is a core for the self-adjoint 
operator . As detailed in [2], is non-negative, and has a one-dimensional kernel spanned by the 
ground state 


2 





In this work, the main question we address is the existence or nonexistence of a non-vanishing spectral 
gap above the ground state. For any Hamiltonian H, such that H > 0 and 0 G spec(i/), the spectral gap 
'y{H) is defined as follows: 

( 6 ) l{H) := sup{(5 > 0 : spec{H) fl (0, S) = 0} 

with the convention that 7 (i?) = 0 if the set on the RHS is empty. In the latter case we call the model 
gapless. We will often denote 7 (i^^) by 7 (A) and 'y{H^) by 7 ( 1 ?) to simplify the notation. 

2.2. Summary of Results. We recall that the domain D depends on a unit vector m G and that 
the Hamiltonian additionally depends on a vector of parameters A = (Ai,..., Xd) G (0, 00 )'^. We will 
often use the notation log A for the vector (log Ai,..., log A^). The goal of this work is to determine for each 
combination of unit inward pointing normal m G and parameters A G (0, oo)"^ whether the gap ^{H^) 
vanishes or not and to find explicit bounds for the gap whenever possible. The case d = 1 was treated in 
[T] where it was shown that model is gapless if A = I and gapped otherwise. Special for the case d = 1 
is that the ground state space for the half-infinite chain [l,oo) C Z is two-dimensional if to log A < 0 and 
one-dimensional if to log A > 0. From now on we will only discuss dimensions d > 2, in which case the ground 
state space for any half-space is one-dimensional. 

For d > 2, the positivity of ^{D) is determined by the angle, d, between the vectors m and — log A, which 
is well defined except in the case A = (1,..., 1). It was already proved in [5] that the model is gapless for 
A = (1,..., I) and is not considered here. 

Our hrst result is an upper bound. For its statement, we define 

(7) c{v) := min{|uj| : Vj ^ 0},v G 

Theorem 2.1 (Upper Bound). For all d > 2, Xi,...,Xd G (0, 00 ), and unit vectors to G such that 
TO • log A < 0, one has the following upper bound: 

( 8 ) 7 (g^)< J |l 0 gA|||sin(d)|, 

c{m)c[X)^ 

where 9 is the angle between the vectors —to and log A. In particular, the gap vanishes if 6 = 0. 

The vanishing of the gap at d = 0 is due to the appearance of extended edge states. The probability 
distribution for the position of the particle in the one-particle ground state is proportional to A^^. It is easy to 
see that, if log A is an outward normal to the boundary, the one-particle ground state assigns approximately 
equal probability for the position of the particle everywhere along the boundary, and the probability decays 
exponentially in the distance from the boundary. 

The following theorem shows that whenever the upper bound given in ® does not vanish, the model does 
in fact have a non-vanishing spectral gap. 

Theorem 2.2 (Existence of a Spectral Gap), //log A yf —1| logA||TO, then ^{H^) > 0. 

This theorem states that if the angle between log A and —to is nonzero, then is gapped. 

2.3. About the Proofs. The results in this work are proved by deriving upper and lower bounds on the 
spectral gap. Theorem l2.1l is proved by a variational calculation in Section [51 For the proof of Theorem l2.21 
we use a well-known relationship between the infinite volume spectral gaps and hnite volume spectral gaps 
as well as the martingale method, all of which we described in detail in Section l3.II 

Our emphasis in proving upper and lower bounds is to establish in all cases whether there is a positive 
spectral gap in the thermodynamic limit or not. We have not attempted to obtain best possible bounds. For 
instance, in the gapless cases we have not established the rate with which the finite volume lowest energy 
excitations vanish in the thermodynamic limit. The gap vanishes at least as fast as 0{L~^), where L is the 
diameter of the support of the excitation. In some cases we can show an upper bound of the form 0{L~^). 
We have not shown that the latter behavior holds in general. 

Since the martingale method proves lower bounds for spectral gaps for finite volume Hamiltonians and 
the PVBS models are translation invariant, the following corollary immediately follows from our proof of 
Theorem 12.21 
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Corollary 2.3. The PVBS model with one species of particle on Z'^ and model parameters Ai,..., Xd S 
( 0 , c») is gapped if there exists at least one j such that Xj ^ 1. 

It was shown in [5] that the PVBS model on Z"^ is gapless if all Xj = 1. Therefore, Corollarv l2.3l comDletes 
the gap classification for the one species PVBS models on Z"^. 

The proof of Theorem l2.2l for a general dimension d is more easily understood by first considering the case 
of d = 2. For this reason, we prove in detail Theorem 12.21 for d = 2 in Section IH and refer to it as necessary 
when proving the general result in Section [5] Furthermore, in addition to the statement of Theorem 12.21 we 
also prove explicit lower bounds for the well-chosen finite volumes and, as a consequence, we obtain lower 
bounds for the gap in the thermodynamic limit. The dependence of these bounds on the parameters m and 
A is somewhat involved, however, which is why we did not include the bounds in the statement of Theorem 
12.21 They can be found in Sections 0] and [5] 

To simplify the proofs, we will permute and reflect coordinates in Z"^ such that the components of the 
inward unit normal of the half-space under consideration are non-negative and the hrst component is non¬ 
zero. That this can be done without loss of generality is easy to see based on the following observations. 

First, replacing mj by —nrij, corresponds to a reflection of A through the hyperplane normal to the 
jth basis vector Cj. The only terms in the Hamiltonian affected by this reflection are the interactions of 
nearest neighbor pairs of the form (x, x + Cj). Let R denote the unitary interchanging the tensor factors in 
Rx ® TLx+ej- Then ^ and it follows that Rh^^^R* = h^^ \ So, we can make all components 

of m non-negative by replacing some of the parameters Xj with X~ ^ . 

Second, since we can assume that all components of m are non-negative and since m 7 ^ 0, at least one 
of the components is strictly positive. We can therefore relabel the coordinates such that mi > 0. Such a 
relabeling corresponds to a permutation of the parameters Xj. 

For future reference we summarize these observations in the following remark. 

Remark 2.4. The model with parameters Ai,..., Ad € (0, 00 ), defined on finite subsets A of the half¬ 
space D C , bounded by the hyperplane containing the origin with inward normal m = (mi,... ,md) G 
is unitarily equivalent to the model on a finite volume A with parameters Ai,..., Ad G (0,oo) and normal 
vector rh such that rhi > 0 , and m 2 > 0 ,..., rhd > 0 . 

3. Iterative Method for Proving Lower Bounds for Spectral Gaps 

3.1. The Martingale Method. To prove that there is a non-vanishing spectral gap above the ground state 
in the thermodynamic limit we rely on the following well-known theorem to reduce the problem to finding 
lower bounds for the gaps for a suitable family of finite volumes. 

Theorem 3.1. Let be the GNS Hamiltonian associated with the connected infinite volume D with spectral 
gap 'y{D). Then for any sequence of increasing and absorbing volumes D, 

j{D) > limsup7(Ai), 

where '){A]fi) is the spectral gap of the frustration-free Hamiltonian H^^. 

To estimate finite volume gaps, we use an approach called the martingale method which is given in the 
following theorem. It provides conditions under which the spectral gap for a frustration-free model on a 
finite volume A, such as the PVBS model considered here, can be bounded by a fraction of the gap for the 
model on small subvolumes. 

Theorem 3.2 (Martingale Method, |10]b For a finite volume A and frustration-free Hamiltonian H^ let 
An be a finite sequence of volumes with Aq = 0 and An A^ = A such that the following three conditions 
hold for the local Hamiltonians for the same i > 2: 

(i) For some positive constant de, 

L 

n—i 

(ii) For some positive constant 7 ^ and nt, if n> Ui, 

ijAn\A„-£ > ^^(l_ 
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where is the orthogonal projection onto Qa = ker(iJ^). 

(Hi) There exists a constant eg < -^ and ni such that ng < n < L — 1 

where En := — G^^+i is the projection onto C/a„ H ^/a„+i ■ 

Then the spectral gap of satisfies 

7(Al) > 

In a typical application to quantum spin chains, the sequence (An) is simply a sequence of increasing 
intervals, say A„ = [1, n], and the conditions are satisfied with dg = £, and £ = 2. Spin ladder systems can be 
treated in the same way with A„ = [1, n] x [1,2] C Z^. In [2] we applied the method to PVBS models defined 
on d-dimensional boxes of the form [1, Lg] x • • • x [1, A^], by using induction on the dimension. For example, in 
the case that d = 2, we apply Theorem l3.2l to the sequence A„ = [1, Li] x [1, n], n = 2,... L 2 - Condition (ii), 
with = 2, is shown to be satisfied by applying the theorem to the spin ladder A„\A „_2 = [1, Ai] x [n— 1, n]. 
Due to the translation invariance of the model this yields a lower bound for the gap for the spin ladders 
independent of n. Applying the martingale method to the spin ladder provides a lower bound of the spectral 
gap independent of Ai and A 2 . This reasoning can be repeated to estimate the gap for d-dimensional boxes. 
The result is a lower bound of the general form 

(9) 7(2'^) > 2-'^7(Ad)(l - V2eif ■■■{!- V2edf, 
where Bd is the d-dimensional unit hypercube. 

In the situation at hand, we will prove a lower bound on the gap for a sequence of finite volumes A^ that 
increases to the half-space A) as A —>■ 00 . Unless m = Cj for some j = 1,... ,d, Aa cannot be rectangular 
boxes, i.e., a cartesian product of intervals. We will choose volumes Aa that are bounded by 2d hyperplanes, 
one of which aligns with the bounding hyperplane of A). In order to avoid a vanishing gap along the 
sequence Aa, the outward normals of all bounding hyperplanes must have a non-vanishing angle with log A. 
This requirement follows from Theorem 12.II The size of Aa will scale linearly with A in each direction. 

For each Aa, we will construct a suitable sequence A„ for which the martingale method applies. In analogy 
to the treatment of rectangular boxes referred to above, we will apply the martingale method inductively as 
the volumes grows in one coordinate direction at a time. Care must be taken so that the Conditions (i)-(iii) 
can be verified for each step. We discuss what this involves for each condition. 

Condition (i) is given by an upper bound on the number of subvolumes A„\A„_^ that contain the support 
of any given interaction term. In this paper, it will always be clear that we can take dg = £. For this reason, 
we will not refer to dg in our estimates. 

Condition (ii) is slightly complicated by the fact that, depending on the normal m, D may not be trans¬ 
lation invariant in sufhciently many directions. It is therefore not possible, in general, to choose sequences 
A„ such that A„ \ A„_^ are isomorphic for all values of n, as was the case for rectangular boxes. As a 
consequence, the quantity jg for Condition (ii) will be the minimum over a finite number of finite-volume 
gaps. 

In order to satisfy Condition (iii), the volumes A„ \ An-g should be connected (in the nearest neighbor 
sense). As demonstrated in Figure O this requires that we consider larger values of i due to the slanted 
boundaries, and this value may be different for each application of the martingale method. It is important, 
however, that we will need only a finite number of different values of i, i.e., the d application of Theorem 
13.21 can be achieved with bounded £ 1 ,... ,id, uniformly in A. This leads to lower bounds of the form 

(10) 7(Aa) > min 7 (A) 

A6V(^1 £3 

where V{£) is a finite collection of volumes for which the length in the j-th coordinate direction of any 
A £ V(£) is proportional to £j. Furthermore, the set V{£) is the same set for the lower bound estimate for 
every volume Aa . 
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In general, the most difficult step in applying the martingale method is finding a sequence of volumes 
which satisfies Condition (iii). However, when applying the martingale method for the PVBS models on a 
sequence of connected lattices such that each subvolume A„\A„_f is also connected, the operator norm in 
Condition (iii) can be exactly calculated and is written in terms of the normalization coefficients for the one 
particle ground states on several volumes. This result is given in the lemma below, and the proof is given in 
the Appendix. 


Lemma 3.3. For each i > 2, and a sequence of increasing finite volumes A„ K, n> 1, such that An and 
An\An-e are connected for all n, the operator norm in Condition (iii) of Theorem \3.Si avvlied to the PVBS 
model (0) is given by 


( 11 ) 


||^A„+i\A„+i_f p ||2 _ C'(An+l-^)C'(A„-|_i\An) 

""" C(A„)C(A„+i\A„+i_,)’ 


where C(A) denotes the normalization coefficient defined in m- 


It is remarkable that the operator norm can be exactly calculated. This lemma will be applied to a variety 
of situations. In each case we will find a simple upper bound for the RHS of (ED, depending on the geometry 
of the volumes (i.e., A„), and 


3.2. Choosing Volumes. From the above discussion it is clear that we need to make judicious choices for 
the sequences of finite volumes appearing the proofs. As we have already discussed, the following three 
properties are necessary for the finite volumes A^ and A„ \ A„_^. First, the boundary of A„ has to include 
an increasingly large subset of the boundary of D (i.e., the hyperplane with inward normal m). Second, the 
orientation of large boundary surfaces should be such that the angle between their outward normal and the 
vector log A remains bounded away from zero. Third, they should be connected. However, the satisfaction 
Condition (iii) is also dependent on choosing volumes that fulfill the following heuristic. 

The general heuristic for selecting a sequence of volumes such that Condition (iii) holds is to choose the 
volumes where the magnitude of the single particle ground state at site x, = exp(2xTog A), is maximized 
at a unique point in each volume when considered as a subset of W^. The maxima occur on the boundary of 
the convex volumes at the points furthest in the log A direction. The uniqueness of the maximum implies that 
the maxima over A„+i \ A„ and A„ \ A„_i are not equal because there is a unique maximum in An+i \ A„_i. 
Due to the exponential form of the 1-particle ground state wave function, the maxima differ approximately 
by an factor not equal to one and consequently the normalization constants for A„+i \ A„ and A„ \ A„_i 
differ by this factor. It follows that the normalization constants in the RHS of (ITT]) are approximately equal 
to a ratio of geometric sums that decays, which is sufficient to satisfy Condition (iii). 

If there is not a unique maximum. Condition (iii) fails. In this case, maxima differ by a factor approxi¬ 
mately equal to one and the RHS of (ITTl) are not geometric sums but sums of ones. It follows that the operator 
norm of Condition (iii) will converge l/v^ as n —> oo and the martingale method will not give a strictly 
positive lower bound. In this case, the set of maximizing points form a line-segment or bounded hyperplane 
spanned by vectors perpendicular to both log A and the normal vectors to the boundary hyperplanes which 
contain a maximum. As the volumes grow as L —>■ oo, these sets of maxima can support extended states 
which in the finite volumes have energy converging to zero. These low energy extended states are an artifact 
of the choice of finite volumes in the sense that they do not correspond to low-energy states in the infinite 
system we are considering. By choosing volumes which have a unique maximum, we avoid this issue. 


4. Existence of a Spectral Gap in Two Dimensions 


Theorem 4.1 (Existence of a Spectral Gap for d = 2). Let D := {x & : m ■ x > 0}. If 

(logAi,logA 2 ) A -|| log All (mi, m 2 ), 


then is gapped. 


We give a detailed proof of applying the martingale method to where D C to motivate the proofs 
in higher dimensions. Thus, we choose a sequence of finite volumes A^ that increase to D such that the 
martingale method proves a nonzero lower bound for the spectral gap of each Hamiltonian which is 
uniform in L, and appeal to Theorem 13.II to show the spectral gap 7 (D) also satisfies the same lower bound. 
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As discussed in Section [5^ when applying the martingale method to PVBS models with slant boundaries, 
it is essential to make good choices for the geometry of the finite volumes. Let 9 be the angle between —m 
and log A. For d = 2 and 9^0, the general heuristic from Section [32] is satisfied for any sequence of volumes 
for which log A is not an outward pointing normal to any boundary. 

Using the symmetries of the PVBS model outlined in Remark 12.41 we prove the gapped statement for 
d = 2 by considering two types of finite volume sequences. In the case that 0 ^ tt, we use parallelograms 
with a pair of boundaries parallel to the boundary of D. We further divide this into two subcases based on 
the inward normal m. Case la is the case where there is a j such that Xj ^ 1 and ^ 0. Case lb is the 
case where such a j does not exist. The choice of paralellograms fails if 0 = tt as then log A is the outward 
pointing normal to the side of the parallelogram opposite to the boundary of D. In this case, we modify the 
parallelogram to a trapezoid. We once again break this into two subcases. Case 2a is the case that both mi 
and m 2 are not equal to zero. Case 2b is the case where either mi or m 2 is zero, but not both as m is a unit 
vector. 

The process of applying the martingale method is the same strategy for all four cases. We first define the 
sequence of finite volumes A D. We then bound the normalization coefficient C'(A) for a general hnite 
volume of the geometry we have chosen, which we use when applying Lemma 13.31 Lemma 13.31 only applies 
the finite volumes are connected as subgraphs of This condition holds for all volumes we choose except 
for in Case la. There, we will determine a sufficient condition for connectedness. In all other cases, we will 
forgo any additional comments on connectedness. 

For each case, we apply the martingale method twice. After the first application (in the direction of X 2 ), 
we will need to estimate the gap for a sequence of quasi-one dimensional systems obtained from the original 
sequence of hnite volumes. This is done by a second application (in the direction of xi) and yields a uniform 
lower bound. The following function / will often appear in the calculation of Ci for Condition (iii): 


finj) 


A2C-i)(1-A2("+i-^))(1-A2) 

(1-A2^)(1-A2-) 


SpeciHcally, there will be some positive constant C for which we hnd estimates of the form 
(12) <C-f{n,£). 


The function f{n, i) is increasing in n. By treating A < 1 and A > 1 separately and taking the limit n —> c» it 
follows that f{n,£) < min(l, which decays exponentially in £. Hence, there exists a minimal 

£ that depends on C such that 

(1-A2)min(l,A2(^-i)) ^ 1 
1 -A 2 ^ ^ £' 

In these situation, Condition (iii) of the martingale method is satished for 

(13) u = VC’min(l, \Zx 21 

for smallest value £. In any situation where we hnd an upper bound of the form (1121) . we will immediately 
use this dehnition for ei. 


Case la: Suppose that (logAi,logA 2 ) E II log A|| (mi, m 2 ) and there exists a j such that Xj E 1 and vtij E 0 . 
Without loss of generality, permute the indices so Ai ^ 1 and mi E 0- 


Volumes: The sequence of volumes A^ yX' D we choose are parallelograms that align with the boundary 
of D. Other parallelograms need to be considered to verify Condition (ii) of the martingale method so we 
dehne a general parallelogram P{bi,b 2 , Li, L 2 ) using a base point, b = (bi, 62 ), and length parameters Li, 
L 2 . Specihcally, 

P( 6 i, 62 , Li, L 2 ) = {x € : 0 < — ■ {x — b) < Li, 0 < X 2 — 1>2 < ^ 2 }- 

mi 

See FigurelD The sequence of volumes Al E' D is given by Al = P{^L, —L, 2L, 2L). 

To apply Lemma |3.31 we need to ensure our volumes are connected. We prove that the condition Li > 
^ + 1 is sufficient to guarantee the connectedness of P( 6 i, 62 , Li, ^ 2 )- Let X 2 = C 2 and X 2 = C 2 + 1 be two 
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X2 



Figure 1. A general parallelogram P(6i, 62,-^1,-^2)- 


subsequent values between 62 and L2 + 62- The parallelogram P(6i, &2, Ti, L2) is connected if and only if 
there is a a;i = ci such that (ci,C2), (ci,C2 + 1 ) S P{bi,b2, Li, L2). To guarantee the existence of such a ci, 
it is sufhcient to show 

( 14 ) X'r(c2) + l<ic 7 (c 2 + l) 

where x(“™(c 2 ) = —^{x 2 — ^ 2 ) + is the lower bound for all xi associated with X 2 = C 2 , and a;“^(c 2 + 1) = 
— ^^{x 2 + 1 — 62 ) + ^1 + Ti is the upper bound for all Xi associated with a ;2 = C 2 + 1. This is summarized 
in Figure[2] Simplifying the inequality in (IT4l) shows this is equivalent to Li > ^ + 1. In particular, is 
connected if 2L > — + 1, and we start the sequence with a sufficiently large L. 



Figure 2 . Diagram for connectedness. 


Normalization Coefficient: For the parallelogram P{bi, 62 , Ti, T 2 ), we always choose 62 , Ti, and L 2 
to be integer valued, while 61 need only be real. The condition 0 < ^ • (x — &) < Li is treated as a bound 
on xi, where the smallest integer value of a;i in the paralellogram for fixed X 2 is 


min 

Xi 


— m2 

mi 


( X2 - ^2) + 61 + r ( x2 ), 
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for which r(x 2 ) € [0,1) is a remainder term that only depends on X 2 , see Figure [TJ For a connected 
parallelogram P(bi, b 2 , Li, L 2 ), the normalization constant is given by 


(15) 


C(P(bi,b2,Li,L2))= 

xGP(bi,b2,Li,L2) 

hj + La-l bi+Li-l-^(x2-b2)+r(x2) 

= E ^ 2 “^^ E 

X2—b2 Xi—bi — ^(x2—b2)+r(x2) 


L 2 1 —T712 ^ 

= (A“ A2)^"=Af("=) E 

ai2=0 xi—0 


Martingale Method in X 2 Direction: Recall that the sequence Al D is defined by Al = 
P{^L, —L,2L,2L). Let N = 2L and define An'^ = P(^L, —A, 2L, n). Notice that An'^ only increases 

f 2 ) 

in the a; 2 -coordinate and that A)^ = A^. 

A ( 2 ) 

We use the martingale method to find a lower bound of the spectral gap of iL w . The martingale method 
allows freedom in choosing any £2 that satisfies the necessary conditions for Theorem 13.21 In general, there 
will be an infinite number of such £2 that satisfy the conditions of the martingale method. However, we 
always choose the lowest possible value of £ 2 - 

As mentioned in 13.21 Condition (i) is satisfied with = £ 2 - For Condition (ii), it suffices to define 7^3 as 


7^2 := min 7 (A[f^\A|f^» ) 

' ^ l2<n<N 


0 ^t> 2 ^ 2 Z/— ^2 


^2 , 


= min 7 P(- 52 , 62 , 2 L,£ 2 ) 


mi 


which is greater than zero since it is a minimum of a finite collection of spectral gaps for finite dimensional 
Hamiltonians. Since this gap is dependent on L, it could happen that 7^3 —?► 0 as A —>■ 00 . However, in the 
second application of the martingale method, we show that for all A, the set of spectral gaps we minimize 
over to define 7^3 share a common nonzero lower bound. We will apply the same reasoning and strategy for 
the remaining cases. 

We apply Lemma [331 as 2A > m 2 /TOi + 1 guarantees connectedness of all necessary volumes. Define 
A 2 = A 2 A)"™^'^"*F Since log A ^ ±|| logAUm, we find that A 2 = A 2 A)""*^^"*^ ^ 1. It is easily seen that for 
A/c C A^, 

C(A„\Afe) = C(A„)-C(Afe). 

Using this and equation (IT5|) we find 


A„||^ = 


\2n\2r'(n) ^n-i2 \2x2 \‘^r{x2) 

^2 ^1 * 2^x2^0 ^2 


(16) 


E n—1 \2x 2 \‘^r{x2) 1 

X2=0'^2 ‘ I^x2=n-e2 + l '^2 ^1 


2x2 \ 2r{x2) 


< max(Ai, A;^ 


A7-')(l-A7+i-^=))(l-Ai) 

(1-A7)(l-Ai") 


The inequality follows from applying the bound min(l,Aj) < ^ max(l,A 2 ), and noticing 


“ niax(Af, A^^ ^). The upper bound in (ITSl) is of the form described in (fT^ . Therefore we define 


e ^3 = max(Af, A^ min(l, ) 


' 1 -Ai 
1 - ■ 


for the smallest value of £2 such that A .11 of the conditions for the martingale method hold for this 


value of £ 2 , so the spectral gap of is bounded below by 


0^62 ^2Z/— £2 


7 (Al) > min 7 ( P(-— 62 , & 2 ,2A, £ 2 ) 


mi 

9 


(1 - e^ 3 v^)^ 


(17) 








Martingale Method in xi Direction: Fix &2 and L. We apply the martingale method in the direction of 
xi to P(—^62, &2, 2 L, £2) to obtain a lower bound on 7^2. Let iV = 2 L and define A,:= P(—^62, &2, n, £2)- 
For these volumes to be connected, we require n > — + 1 . For the same reason — + 1 also serves as a 
lower bound for £1. 

Condition (ii) is satisfied with 7£j^ defined by 


:= min 

li<n<N Ti *1/ 


min 

-^b2<&i<-L-4-^b2 


7(P(6i,62,£i,£2)) 


For any choice of (61,62), the area of the parallelogram P(6i, 62,-^i, £2) is ^1^2- If the li are fixed, and 
•^1 > ^ + 1, then there are only a finite number of distinct connected subsets of 1 ? (up to translations) 

that are contained in a parallelogram of area •^i-^2- Let P(£) denote the collection of all such distinct subsets. 
Since the PVBS models are translation invariant, it follows that 


( 18 ) 


min 7(P). 
Pevie) 


By definition, ^{i) is nonzero and independent of 6 and L. Therefore, we can use 7)^) for Condition (ii) in 
place of 7£j. 

For Condition (iii) recall that, by assumption, Ai ^ 1 . Using similar calculations and bounds of the 
operator norm as computed in the X2 direction, we find 




\ 2 n ,2x1 

^1 Z^xi=0^1 


E ra—i Y^ra—i ,2 

xi=0^1 ■ Z^xi=ra-£i + l 

^ • /, ,2(^i-l), 1 — ^1 

<min(l, Ai' ')-- 

i - Ai 


, 2 X 2 


Which is also of the form of (TT^ . Therefore, define 
( 19 ) := min(l, Ai^“^), 


1 -A? 


1 - Af 1 


for the smallest value of £1 > ^ + 1 such that Then all conditions for the martingale method are 

satisfied for this value of £1 and 7(A^^) satisfies the nonzero lower bound; 

(1 - 

i\p, 

PeP{e) 


( 20 ) 


7(Aw^) > min^7(P)i 


Since this result holds for all choices of 62 and L, it follows that this is also a lower bound for 7^^, which 
produces a uniform lower bound on 7(Ai). Applying Theorem 13 . li on 7(P), we find 

7(P) > mm 7(P)^- - -^ > 0. 

P^V{i) *2 


Case lb: We still consider (logAi, logA2) ^ || logA||(mi, m2), but now we assume that for all j either 
Xj = 1 or rrij = 0 . Since we are also assuming that (logAi, logA2) ( 1 , 1 ), as this is a gapless case, and 
(mi, m2) ^ (0, 0), as m is a unit vector, this falls into two cases: 

Ai 1, A2 = 1, mi = 0, m2 = 1 or Ai = 1, X2 ^ 1, mi = 1, m2 = 0. 

We motivate the need to consider these cases separately. In the former case, the definition for the slant 
boundaries of the parallelograms P(6i, 62, Pi, P2) are undefined. In the latter case, the parallelograms 
P(6i, 62, Pi, P2) are rectangles. Since in this case, (log Ai, log A2) = c( 0 , 1 ) is the outward pointing normal 
to some boundary of P(6i, 62, Pi, P2), this is not a good sequence of volumes to apply the martingale 
method. 
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Volumes: By Remark 12.41 it is sufficient to only consider the case Ai 1, A 2 = 1, mi = 0, m 2 = 1. For 
this choice of m, the infinite volume D is the upper half plane {x € : X 2 > 0}. Since (logAi, logA 2 ) = 

c(l, 0), we must choose finite volumes with no vertical boundaries. So, we consider parallelograms with the 
vector V = (1,1) generating the slant boundaries. Specifically, 

P{bi,b2, Li,L2) = {x € I? 0 < V ■ {x - b) < Li, 0 < X2 - &2 < L 2 } 

We choose = P{—L, 0, 2L, 2L) to be the sequence of volumes increasing to D. This is the simplest 
sequence of parallelograms to choose as given bi, & 2 : Ai, L 2 are integer valued, the smallest and largest 
values of Xi value for any X 2 will lie on the boundary. Hence, there are no remainders r{x 2 )- As such, the 
proof for this case is the similar to Case la, with several simplifications by our choice of v. 

Normalization Coefficient: Since we no longer have a remainder r(x 2 ) and since A 2 = 1, the normal¬ 
ization constant calculated in (1151) can be simplified to 

C{P{bi,b2,Li,L2)) = A^" £ £ Afb 

X 2 —O xi—0 


Martingale Method in X 2 Direction: Define N = 2L and let h.n'^ = P{—L, 0, 2L, n). Then = A^^^ 
For any value of £ 2 , our choice of v is such that the volumes An\An-e 2 are isomorphic sublattices of Z^. By 
the translation invariance of the PVBS models, Condition (ii) of the martingale method is satisfied with 

It only remains to check Condition (iii). Using Lemma 13.31 summing the geometric series and simplifying 
the expression yields 

" " “ (i_Ar“)(i-Ar"=) 

which is of the form m- Therefore, we define 


ee, = min(l, A^ 


It can be shown that the smallest value of £2 such that is £2 = 2. Choosing this value of £2 gives 

min(l, A^^) 


£2 = 


1 -l- Aj^ 


-2 


Thus, we have the following lower bound on 7 (Ai) : 

( 21 ) 


7(Aa) > A(Af). (Ljaklt 


where = P{-L, 0, 2L, 2). 

Martingale Method in xi Direction: We now apply the martingale method to A^^^ to obtain a lower 
bound on 7 (A 2 ^^) that is independent of L. Let N = 2L and define A^^^ = P(—L,0,n,2). Similar to the 
appliation in the X 2 direction. Condition (ii) of the martingale method is satisfied for 

7 ^, =7(Alf) = 7(^(0,0,£i,2)). 

For Condition (iii), a similar calculation as the application in the direction of X 2 shows 

(1 - Af)min(l, 


(22) 




1-Af 


2^1 
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This is the same bound we had for the first application with replaced with Ai. Therefore, Condition 
(iii) once again holds for and 

min(l, Ai) 

\/l + ^1 

Thus, we see that 

(23) ^(A(2))>^(P(0,0,2,2))^i^l^^. 

Substituting this lower bound into (I^Tl) produces a lower bound on 7 (Al) that is independent of L. 
Therefore, we obtain 

(24) ^(i7)>7(P(0,0,2,2))-li^l^^.ii^^^ >0. 

Case 2a: We now consider (log Ai, log A 2 ) = || log A|| (mi, m 2 ). We invoke the reflection symmetry of the 
model to produce a boundary for which mi > 0 and m 2 > 0, see Remark 12.41 For this case, we also assume 
that mi ^ 0 and m 2 ^ 0. Since log A = || logA|| • m, this implies that Ai > 1 and A 2 > 1. 

Volumes: We cannot choose parallelograms for the sequence Aji, since (log Ai, log A 2 ) = c(mi, m 2 ) would 
be an outward pointing normal to the boundary opposite to D. We instead use trapezoids which have one 
boundary along the boundary of Z?, and vertical or horizontal lines for the other boundaries. The general 
trapezoid of interest is defined by 

T(6i, 62, ^1,^2) = {a; e : 0 < — • {x - b), xi - bi < Li, Q <X 2 -b 2 < L2) 

mi 

where b = (&i, 62 ) is the base point of the trapezoid and Li, L 2 are the lengths in the coordinate directions, 



Figure 3 . A general trapezoid T(6i, 62, ^1,^2). 

see Figure[3] We always choose 62 , Li, and L 2 to be integer valued. However, bi will typically not be integer 
valued. For the sequence Al D, we choose Al = T{^L, —L, 2L, 2L). 

Normalization Coefficient: For a fixed integer value of X 2 , the smallest integer value x™™ such that 
(x("™, X 2 ) is contained in the trapezoid r( 6 i, 62 , Ai, L 2 ) is given by 

xr" = - — {X2 - 62) +bi+ r(x2) 
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where r{x 2 ) G [0,1). In particular, since bi need not be integer valued, it follows that a;™™( 62 ) = bi +r{b 2 )- 
For all X 2 , the largest integer value of xi, denoted contained in the trapezoid T{bi, 62 , Li, L 2 ) is 

x-“" = Li + 6 i-l + r( 62 ). 

Since logAi = cmi > 1 and log A 2 = cm 2 with c > 0, it follows that A 2 • ^ 


^2 +-£^ 2 ” 1 -L 1 +&I — l+r(& 2 ) 

C(T(6i, &2, ii, i2)) = ^ A 2 "= ^ 

2 : 2=62 xi=bx-'^{x 2 -b 2 )+r{x 2 ) 




L2-I . 2(Li+i’(62) + ^a:2-i'(x2+b2)) 

^ y 2 b ^ ^2r(a;2+62) , K _ ZJl 


X2—O 


Af.-l 


C - ^ and C ■ - 

^max — x2 1 cxiiu. ^min — ^2 1 ' 

— i — i 


(25) 

Define 

We show that 

,'^i^l+r{b2)+^X2-r{x2 + b2)) 

/r,a\ r< ^2(Li+r(62) + :^a;2-r(a:2+62)) ^ A^ ~ ^ ^ r< ^ 2(Li+r(b2) +:^X2-r-(a:2+62)) 

(zbj S: ^^2 

The upper bound follows immediately. For the lower bound, note that 

Li + r{b2) + —X2 - r{x2 + 62 ) = Ti + <“(^ 2 ) - 

mi 

where a;™*" is the minimum xi value for X 2 + 62 - Since toi, m 2 > 0, as X 2 increases the minimum value of 
xi decreases. Therefore, a;™™( 62 ) — a;™*" > 0. Trivially Li > 1, and so 

1 < Ti + r{b2) H-^a ;2 — r{x2 + 62 ) 

mi 


and the lower bound follows. Inserting these bounds into (l25l) yields 


L 2 —1 


X/ 2 ~ 1 


(27) ^ A^"^ < C{T{bi,b2,Li,L2) < ^ A^"^ 


312=0 


ai 2=0 


Martingale Method in X2 Direction: Let N = 2 L and set An'^ = T{^L, —L, 2 L, n). Then for n > £2, 

= T(- —(n - L - £2), n - L - £2, 2L + —(n - £2), £2). 

TirL\ T71j\ 

Therefore, Condition (ii) of the martingale method is satisfied for 

7^3:= min 7 ( r(- —(n - T - £2), n - L - £2, 2L + —(n - £2), £2) 

l.2<n<N \ mi mi 

= , 7 (?"(- —&2, &2, 2 L+ —(L + 62), £2) 

-L<b2<2L-l2 \ mi mi 

For Condition (iii), notice that for all fc < n 

Ai2)\A^2) ^ T{- — {k -L),k- L, 2 L + —k, n - k). 
m\ mi 

Using (EZl) for each of the four normalization constants we compute 

- 7 “- (i-Ar=)(i-Ar)— 

This is of the same form as m- Since A2 > 1 , we choose 


(29) 


Q, := 


Cmax / 1 -^2 


(7 . V 1 ’ 

^min y -L — ^>2 
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for the smallest value of £2 satisfying < j^. Therefore, Condition (iii) of the martingale method is also 
satisfied and 


(30) 


"/(Al) > min 

-L<b2<2L-e 


T{ - -b 2 , & 2 , 2L H- {L + 62 ), ^2 

TOi mi 


(1 - 

(2 


Martingale Method in xi Direction: We apply the martingale method to produce a lower bound on 
7 ^ 2 . Fix L and 62 , let TV = 2L + + 62 ), and define = T(-^ 62 , ^> 2 , n, £ 2 )- 

For Condition (ii), recall that the martingale method assumes that = 0. So, for any £i > 2, 

AeMo=T{-—b2,b2,ii,i2). 

mi 


7 t(/) = min 7 (r( 6 i,& 2 , •^1,^2))- 


Let 

(31) 

Similar to the parallelogram case, there are only a finite number of distinct (up to translations) trapezoids 
T ( bi , b2 ,£ i , i2 ) contained in Therefore, by the translation invariance of the PVBS model, the minimum 
in m is positive, and 7 (A£j\Ao) > 7 t(^) > 0. Furthermore, for n > 1 

An+^i\An = B{ -62 + n — £1,62, f’l, -^2) 

mi 

where i?( 6 i, 62 , ^ 1 ,-^ 2 ) is the rectangular box 

B ( bi , b2 , ii ,£ 2 ) = {x € I ? ■. bi < X < bi + b2 < X2 < b2 + £ 2 )- 

Regardless of the values of &i and & 2 , T3(6i, 62 , ^ 1 , ^ 2 ) — T?(0,0, £ 1 , £ 2 )- Therefore, 7 (A„+^^\A„) = 7 (^( 0 ,0,£i, 
and Condition (ii) of the martingale method is satsified for 

(32) 7 ^j = min (^7 t(£), 7(R(0, 0, £ 1 , £2 
For Condition (iii), since for all n > k, A„\Afc is a rectangular box 

C(AWAAi'^) _ A?('^-')( 1 -Af) 


(33) ,,, ,,, 

c(aWaaWi_,J 

Using (Eli), we additionally obtain the bound 


1-A 


2(1 


(^(Ai^^) Cmin 

C 1 _ X 2 

Aymax ^ /'I 


(34) 

Putting these together yields 

(35) 




£7 1 ^ 2^1 ' 

'^min i — A-^ 


which decays exponentially. 

Choosing £1 to be the smallest integer such that satisfies all the conditions for the martingale 

method and 

(1 - e^iAA)^ 


(36) 7 ( T{-—b 2 , & 2 , 2L + —(L + 62 ), £ 2 ) ) > niin 

mi mi 


,7(i?(0,0,£i,£2 


£1 


> 0 . 


Since this bound is independent of the choice L and 62 , replacing this for 7 ^^ in the lower bound for 7 (Ai) 
and applying Theorem proves 13.11 

(1 — QaV^)^ (1 — Qia/A)^ 


'y{D) > min 


, 7(i?(0,0,£i,£2 


> 0 . 


Case 2b: We now consider (log Ai,logA 2 ) = || log A||(mi,m 2 ) and either m = (1,0) or m = (0,1). These 
cases are considered separately from the those covered in Case 2a for the following reasons. In the case that 
m = (0,1), the trapezoids r( 6 i, 62 , Ti, L 2 ) are undefined. If m = (1, 0), the boundary of D becomes a vertical 


£2)) 
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line, and the trapezoid r(&i, 62 ,-^ 2 ) is the rectangular box i3(6i, 62 ,-^ 1 , ^ 2 )- Since logA = ||logA||TO is 
an outward pointing normal to one of the boundaries of & 2 , ^i, ^ 2 ), this is not a good choice for the 
volume sequence for the PVBS model. 

Permuting the coordinate indices, we need only prove this special case m = (1,0). This implies Ai > 1 
and A 2 = 1 . 


Volumes: We modify the trapezoid T{bi,b 2 , Li, L 2 ) from Case 2a so that the right most boundary is 
slanted with slope determined by the vector v = (1,-1). Specifically, the trapezoids are defined as 

( 37 ) T(6i, 62, Li,L2) = {x £ : 0 < xi — bi, v ■ {x — b) < Li, 0 < X2 — b2 < L2}. 

Since the boundary of D is defined by xi = 0, we choose as our sequence of increasing volumes Al = 
r(0, —L, 2L, 2L). In this case, we only consider trapezoids for which 61 , 62 , Ti, and L 2 are integer valued. 


Normalization Coefficient: We bound the normalization constant for the one particle ground state on 
T( 6 i, b 2 ,Li, L 2 ). Since all 61 , & 2 , Li, and L 2 are integer valued, there are no remainders. For a fixed value of 
X 2 , the smallest value of Xi is always x™™ = 61 , and the largest value of xi is = 61 + + X 2 — 62 — 1 - 

Using A 2 = 1, we obtain 


^2 — 1 .. _ \ 2 (Li+a: 2 ) 

C'(r( 6 i, 62 ,Ui,i 2 )) = A^^ ^ 1 


312=0 


1-A? 


Since Ai > 1, similar to Case 2a, we have 

L2-1 

(38) C™„A2^Af ^ ^ A?"= < C(r( 6 i, 62 ,Ui,U 2 )) < C^a.A^^Af ^ ^ A 


Lo — l 


2X2 
1 ’ 


ai 2=0 


ai 2=0 


where Cmin and Cmax are dehned as they were in Case 2a. 


Martingale Method in X 2 Direction: Let N = 2L and set = T(0, —L,2L,n). Then An^ ■ 

For k < n, 

(39) ^ k-L,k + 2L,n-k). 

Therefore, Condition (ii) is satisfied with 

7^2 = min '){A^^\A^^}_n) 
e 2 <n< 2 L \ n t 2 ' 

= min 'y{T{0,n — £2 — L,n — £2 + 2L,i2)) 

^2<71'<2L 

= min 7 (r( 0,62 - L, 62 + 2 L,£ 2 )) 

0<b2<2L-l2 


,(2) ^ A (2) 


Using 


we find 


O ) ^max ^ ^max 2(.^2 —1) _ 2^2 

^ Ai — to 


Af"- 1 


Cr, 




c{h^:U-,2) 

Since the one particle ground state normalization coefficient increases as the lattice increases we know that 

C(Ai^jAAl^^) 


(40) 

Therefore 


C(Ai^|AAi^|i_,,) 


< 1 . 


|| GA '^ iAAA .-..£;„||2 < ^ aA ^= -.= 4 . 


(41) 
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which decays exponentially. Letting £2 be the smallest integer such that the martingale method 

yields 

( 42 ) 7(A,.)> min 7 (r(0, 62 - L, 62 + 2 L, £2)) • ~ ■ 

0<02<2L—^2 


Martingale Method in xi Direction: We apply the martingale method to the finite volume 
T( 0 , 62 — L, 62 + 2 L, £2) for a fixed L and 62 to obtain a lower bound on 7^3. 

Let = 62 + 2 L, and defined = T{b2 + 2 L — n, &2 — L,n,£2)- Then A^^ = T( 0 , 62 — A, 62 + 2L,£2) 

as desired. Similar to Case 2 a, since Ag^^ is the empty set, we have 

aW\a|,^) = T(&2 + 2 A - £1, &2 - a, £1, £2) = A( 0 , 0 , £1, £2), 
where the last equivalence holds since 62 and A are integers. For n > 1 

Ai\\AW = A(&2 + 2 A - n, 62 - A, £1, £2) ^ A( 0 , 0 , £1, £2). 

Using the translation invariance of the PVBS model. Condition (ii) is satisfied for 
(43) 7^1 = min (7(r(0, 0 , £1, £2)), 7(^(0, 0 , £1, £2))). 


For Condition (iii), using A2 = 1 and C'(A„\Am) = C'(A„) — C{Ak) for n > fc > 1 , we can exactly calculate 
that 

& 2 + 2 Z/—fe—1 - \2{n—k) 

C(A<»\Al,‘>) =1:1: >'"' = 

X 2 —b 2 —L xi—b2-\-2L—n 

It follows that, 

1 - A? 


1 -Af 


c(aWaa1 '^) 


2^1 ■ 


C(AWAAi'|i_,,) l-AA 

Since C Ai^\ it follows that C'(a£^^_^ J\C'(Ai^^) < 1 , so 

\ n 1 — A? 


2^1 


1 - A^ 

which decays exponentially. Let £1 be the smallest integer such that 

2 1 - A? 1 

el :=-< —. 


2^1 


(1 - ei^y/l[y 


1 -A. 

Therefore the three conditions of the martingale method are satisfied and 

7(r(0, &2 - A, 62 + 2 A, £2)) > min (7(T(0, 0 , £1, £2)), 7(-B(0, 0 , £1, £2))) 

Since the choice of £1 is independent of A2 and 62, substituting this into the bound for 7(Ai) gives 
( 44 ) 7(AA > min (7(r(0, 0 , £1, £2)), 7(^(0, 0 , £1, £2))) 


£1 

Since Al Z' D, the same bound holds for the spectral gap of H^, as desired. 


5 . Existence of a Spectral Gap in d Dimensions 

In analogy with the analysis for d = 2 in the previous section, we divide the proofs for arbitrary d into 
several cases. Cases 3 a and 3 b cover the situations in which the vectors m and log A are not parallel, that 
is log A 7^ ±11 logAllm. For these cases, we choose volumes with d pairs of parallel boundaries, i.e. the d 
dimensional analogue of parallelograms, sometimes called parallelotopes. Case 3 a is the case where there 
exists j such that rrij ^ 0 and Aj ^ 1 and Case 3 b is the case where there does not exist such j. Case 4 
deals with log A = || logAHm. For this case, we choose the d dimensional analogue of the trapezoids in Case 
2 a and Case 2 b. We will refer back to arguments and calculations in the proofs for two dimensions. 

The proof of Corollarv l 2 . 3 l is given at the end of the section. 
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Case 3a: Suppose log A ^ ±|| log X\\m and there exists j G {1,..., d} such that A^ ^ 1 and rrij ^ 0. By 
Remark 12.41 we can permute the indices so Ai 1 and mi ^ 0, without loss of generality. 


Volumes: As in Case la, we choose one pair of boundaries given by 0 < ^ • a; < L. If we chose 
—L < Xk < L for the other boundaries, the volumes may not satisfy the heuristic that A^ is maximized 
at a single point. This occurs exactly when (log Ai, log Aj) = c{mi,mj) as the ground state coefficients 
A“ are maximized along some (xi,Xj) plane contained in the volume. To avoid this issue, we introduce a 
different pair of slanted boundaries. The non-parallel condition ±|| logA||m implies there must be Xj such 
that Xi ^ ^ argument is as follows. 

The collection of vectors {— 771 ^ 6 ^ -|- miei : j = 2,..., d} is a linearly independent set such that all vectors 
are orthogonal to m. Consequently, this set spans the perpendicular subspace of m. If log A is not parallel to 
m, then it must have a non-zero projection to this perpendicular space. Therefore, the inner product of log A 
and —mjCj+miei is non-zero for some j = 2,..., d. Then Xi ^Qg jQg ^ i 

We permute the indices so that this holds for j = 2. 

The extra pair of slanted boundaries is —L < v ■ x < L where the vector u = (0,1, U 3 ,... is defined as 

[ 0 if a7^a,- ^ I 

Vj := < " 

i -I if Ai A, = I 

We apply the martingale method to the sequence of volumes is given by 

Al := {a: G Z'^ : 0 < - - x < L, —L < v ■ x < L, and for j > 3, —L < a;, < LI. 

TOl 

Normalization Coefficient: To abbreviate notation, define 

Ai:=Ai, A2:=A7^A2, and Aj := (A”^Aj)(A”^A 2 )“’'^ for j > 3. 


By the choice of Vj, we guarantee that each Aj ^ I for j = I,..., d. We calculate the normalization constant 
for general volume 

TTl 

A' := {x € D : bi < — • x < 61 -|- Li, b 2 < v ■ x < b 2 + L 2 , for j > 3, bj < xj < bj + Lj}, 

which includes all possible subsets generated by the martingale method. For an integer choice of 62 , the 
compound inequality 62 < x • x < &2 + L 2 provide integer bounds on the X 2 coordinates of the form 

d d 

62 + '^{-Vj)Xj < X 2 < 62 + dv 2 - I + 

3=3 3=3 

The compound inequality bi < -^ ■ x < bi + Li bounds the xi coordinates. Equivalent integer bounds are 

given by adding a remainder term r such that ^ equals the integer Note 

that r depends on m and X 2 ,..., x^, but not xi. Hence, the xi values are bounded by 

d d 

61 -I- >- -Xj + r < Xi < bi + Li + y - -Xj + r. 


To simplify notation in the calculation of the normalization coefficients, let ai(x) = 

0 - 2 (x) = X]j =2 ~^^3 1“ check that the normalization coefficient for A' is then given by 


(45) 


C(A0 





yi=i j xd=o xd_i=o x2=o xi=o 

This is bounded above and below by a constant multiple of a product of geometric series after bounding the 
remainder term appropriately. 
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Martingale Method in Each Direction: We apply the martingale method d times. We hrst apply it 
to the sequence of volumes 

■= Al n {xd < U - L}, 

for id < n < 2L to generate a lower bound of the form 

7 (AW) > min 7 (aW \ 

-L+id<n<L “ td 

Introducing bd = n — id, we apply the martingale method to each \ a[^^ using the sequence 

A(-^-i)( 6 d) := \ aIJ) n {xd-i <n-L} 

= Al n {hd <Xd<hd + id} n {xd-i <n - L}, 


which bounds the gap of each \ MfJ by 


Ad) 




(d-l) ,, (1 - ^d-l\/id-lY 


i -1 


We iterate this process d times: the sequence for j-th application martingale method is 
AW , 6d) := \ n {x, < n - L} 

= Al n {for k > j, hk<Xk<bk+ 4} n {xj <n- L}, 
for ij < n < 2L. The last two sequences differ slightly. The sequence in the X 2 direction is given by 

.■.,bd):= \ A[^^) n {u • X < n - L} 

= Al n (for k > 2, bk < Xfc < 4 + 4} fl (u • x < n — A}, 


for i 2 < n < 2L. For the x\ direction, we use the sequence 

A^i)(& 2, ..., 4) := (Aif+^^ \ A[^)) n • X < n} 

Tfi 

= Al n (for k> l,bk<Xk <bk+ 4} n {-x < n}, 

mi 

for ii < n < L. 

Condition (ii) is satisfied by letting := min„ 7 (A^^ \ To apply the Lemma [3.31 for Condition 

(iii), the volumes A^^^ \ A^^2f. intersected with the lattice must be connected. We verify this by checking 
that cross-sections of the volume are graph connected in each coordinate direction. The cross-sections of the 
volume parallel a general (xj, Xk) plane with j,k ^ 1,2 are rectangles and connected. The cross-sections with 
the X 2 coordinate and not the xi coordinate are either rectangles if Vj = 0 or parallelograms with boundaries 
slope one which are connected since L 2 > i 2 > 2. For cross-sections with the xi coordinate, the argument 
in Case la applies under the condition Ai > -^i > -I- 1. We impose the condition L> ii> maxj(^) -|- 1 

to guarantees connectedness of the volumes. 

We now calculate ee .. For the application in the Xj direction, all terms independent of Xj cancel in the 
ratio of normalization coefficients. We find 
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A (3) \aO) 

'■^'•n + 1 V^n + l-i 


En\\ < niax(A{, 


-4 


)• 


^2{e, 

(1-Af )(1-Af) 


which is of the form stated in (jl2F Therefore, we define 


ti- := max(Ai, A;^ min(l, A^^ 


\ 


1 - Xf 
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for the smallest value of such that e^. < j-. By the martingale method, we have 


(^i+1: min7(AW(6j+i,. ..M)\ ■ ■ • > ^d)) 


(f) 




Combining the Martingale Methods: We have defined the sequences such that {bj+i, ■ • ■, bd) = 
\ • ■ •) ^d)- We iterate the method for each coordinate direction by 

applying the next application to a general '){Kn\bj+i,... ,bd) \ ... ,bd)). Each application 

of the method results in a new constant factor and a minimum over sets with fewer degrees of free¬ 
dom. In the final application of the martingale method, Condition (ii) requires a minimum over all 
■ ,bd) \ ..., bd)), which are d-dimensional volumes with coordinate lengths indepen¬ 

dent of the sequence. Let £ denote the vector {£i,... ,id)., and define 

Ar := {x G : 0 < ■ x < £i, 0 < v ■ x < £ 2 , for j > 3, 0 < x, < £A. 

^ nil 

Then, the finite volumes for which we need to lower bound the gap, after a suitable translation in Z'^, are 
subsets of A^. We call the resulting set of translates V{£). The family of sets V{£) itself is finite up to 
translations. Therefore, combining bounds we obtain 

7 (Al) > min^ 7 (P) ^ q_ 

PG-P(f) £j 

Case 3b: Suppose log A 7 ^ ±|| logA||m and that there does not exist k such that Afe 7 ^ 1 and nik 7 ^ 0. 

Recall that we are only considering X ^ (1,1,...,!). Therefore, there is at least one k such that Afe 7 ^ 1. 
Since m is a unit vector, not all mi,..., rud are equal to zero. This implies there must be at least one Afc = 1. 
We permute the variables so that Ai = ... = Xj>-i = 1, and Xji, Aj'+i ,... ,Xd are not equal to one. We 
assume, without loss of generality, that mi = maxj{mj}. 

The proof for this case is essentially the same as the proof of Case 3a. We indicate how to choose the 
volumes A^, how to obtain and £ 1 ,... ,£d, and give adapted definitions for Xj, A^and 'P{£). 

Volumes: The finite volumes for this case are defined as 


Ai = {x G Z'^ :0 < — ■ X < L, —L < v ■ x < L, for j = 2,..., d — 1, —L < x,- < L}, 

mi 

where v = (—1, —2,..., —2,0,..., 0, mi), where Vk = —2 for fc = 2,..., j' — 1. This vector determines two 
slanted boundaries which guarantee the volumes satisfy the heuristic of Section [ 

Normalization Coefficient: Define 


Xj 


Xd 


A 

Ai 


2—mj /mi 
d 


if j = 1 
if 2 < 7 < 7 ' — 1 
if f <j<d 


Each modified parameter satisfies Xj 7 ^ 1. To see this, note that by our original choice of ordering, Xj 7 ^ 1 
for j = j',...,d. Consequently, it also follows that Ai = Ad 7 ^ 1. Since mi = maxj{mj}, it follows that 
nij/mi < 1, and therefore Xj = A^ ™j/™i ^ p 

The general volume for all d applications of the martingale method is of the form 


A' = {x G Z*^ :bi < -• X < foi -I- Li, bd < v ■ x < bd + Ld for j = 2,..., d — 1, bj < Xj < bj -b Lj}. 

nil 

The compound inequality 61 < ^ • x < 61 -I- Li once again bounds the values of xi, where as the compound 
inequality bd < v ■ x < bd + Ld bounds the values of Xd- To abbreviate notation, let ai(x) = J^jS^ “^ 7 ^ + ^ 
and ad{x) = The remainder term for ai(x) only depends on X 2 ,..., x^'-i. The normalization 
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coefficient is then 


C(A') 


bd—i-\-Ld — i — l 62 +-£/ 2~1 fci+ Li — 1 +ai (tc) bd-\-Ld — ^-\-cid{x) 

E E E E 

Xd-i—bd-i X2—b2 xi—bi-\-ai{x) Xd—bd-\-ad{x) 



Ld — ^ 

z. 


Ld-i—l 


\ “^Xd — l 

M -1 


Xd^O 


Xd-i^O 


X/2 —1 

' E ^2 

2 : 2=0 


Li -1 


2a:2 ^2r E 


a:i =0 


Martingale Method: For brevity, we define the and ij and refer to Case 3a for the bound. We 
define 


Q,- := max(A^,A^^)min(l, Xf 

for the smallest value of £j such that e^. < j-. We may drop the error term max(A^, A^^) for all j > j' since 
rrij = 0 implies r does not depend Xj. Analogous to previous case, we define A^ as 

A/ := {a; G : 0 < — • a; < ti, Q < v ■ x < Id, for 2 < j > d — 1 , 0 < a:,- < £,}. 

mi 

Once again, the finite volumes for which we need to lower bound the spectral gap, after a suitable translation 
in are subsets of A^. Let Vi^i) denote the resulting finite set of translates. After d iterations of the 
martingale method, as in Case 3a, we have a positive lower bound on the gap of the form 


eel 

2 L' 
j 


1-A 


7(Al) > 


d 

min 7 (P) TT 

Peril) fj. 


(1 - 


We note that the parameters ij, Cj, and the set 'P{t) are different from those in Case 3a. Once again using 
Theorem 13.II we conclude that 'y{D) > 0. 

Case 4: Suppose log A = || log A||m. 

As described in Remark 12.41 we apply coordinate reflections so each rrij > 0 and permute the indices so 
that mi > 0. We proceed by stating how to choose the volumes Al, how to obtain good values for ^ 1 ,..., id, 
and give adapted definitions for Aj, j = 1 ,...,, d, and two analogues of A^ and Vil). 

Volumes: We define a vector v := {1,V2, ■ ■ ■ ,Vd) with 


For the sequence A/, we take 


r — 1 : rrij = 0 

\ 0 -mj ^0 


Al := 


{x€Z‘^: 


m • X 
mi 


> 0, V ' X < 


d 


^+E( 

J=2 


mi 


Vj)L, for j = 2,..., d, —L < xj < L}. 


The summation term on the bound of u • a; is used to guarantee a length of at least L in the Xi direction for 
any X 2 ,... ,a:d in [-L, L). 

Normalization Coefficients for X 2 , - ■ ■ ,Xd directions: To simplify notation, we define Xj as 

Al := Al, and Xj := A™^ ^ for j > 2, 


which are all strictly greater than one since log Aj = || logA||mj and rrij > 0. The general volume A' defined 
below covers the volumes generated by the martingale method for the X 2 , ■ ■ ■ ,Xd directions. 


A' 


{x G Z"* : 0 < 


m ’ X 
mi ’ 


V • X 


d 


<Ai+E( 

i=2 


mi 


Vj)L, for j = 2,..., d, bj < Xj < bj + Lj}, 
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where bj > —L and Lj < 2L. Define 

ai(a;) = XI and bi{x) ='^{^L - Vj{xj + L)). 

j=2 j=2 

Furthermore, let r{x) be the remainder associated with ai(a;), that is r{x) = |"ai(a:)] — ai(a:), and define 
r{L) be the remainder term associated with X 2 = X 3 = ■ ■ ■ = Xd = —L. The normalization constant is given 

by 


^d+^d —1 

C'(A')= X 

Xd=bd 


b2+L2 — l Li — l+bi{x)+r(L) 

■ X 

X2—b2 x\—ai{x)-\-r{x) 


b2+L2 — l _T, 


^d+^d —1 -'^d - - 

= Y. X (A~A2f"=Af 

Xd=bd 2 : 2=62 

bd+Ld — i b2+L2 — l ^ 2 {Li+bi{x) — ai{x)+r{L)—r{x)) 

^ ^ ^ ^2r-(2:)Ai_ " 

Xd = bd X2=b2 


Li — l+bi(x) — ai{x)+r(L) — r{x) 


X 

tci =0 


A?-l 


In the above computation, we use that Xj = 1 as log A = || logA||m. Just as in the proof of Case 2a, 
the exponent of Ai is greater than 2 since bi{x) — ai{x) > 0 , r(x) < 1 , and we require that Ti > > 2 . 

Therefore, the bounds from Case 2a hold here, that is 


^ 2 (Li+bi{x) — ai(x)+r(L)—r{x)) 
-^1 



- 1 


^ ^ ^ 2 {Li+bi(x)-ai{x)+r(L)-r{x)) 


^ 2 (Li+bi(x)-ai{x)+r(L)-r{x)) _ 


Af -1 


<C„ 


^ 2 (Li+bi{x)-ai{x)+r{L)-r{x)) 


Therefore, C'(A') is bounded above and below by 




n 




^ 2 {bj-\-L) 


Lj 1 

X^ 

Xi—0 


< C{A') < CmaxX 


2 {Li+r{L)) 


n 

i =2 


/ x 2 {bj-\-L) 

I ' 



Martingale Method in X 2 , - ■ ■ ,Xd Directions: The first d — 1 applications of the martingale method 
follow the same procedure, which we now detail. The last application will be slightly different and will need 
to be considered separately. We apply the martingale method first to the sequence 

A^^^ := Al n {xd < n — L}, 0 < n < 2L. 


Replacing n — L hy bd + id, apply the martingale method to each \ A^^^ using the sequence 

^n~^Hbd) = \ Kf) {Xd-i < -A + n}, 0 < n < 2L. 

We iterate this procedure for j > 2, applying the martingale method to the sequence \ 

using the sequence 


AW( 6 ,+i, ...,bd) = ...,bd)\ AiX'^(^.+ 2 ,..., bd)) n {x, <-L + n} 

on (ii) is satisfied by taking the mir 
may apply Lemma 13.31 and compute the norm 


Condition (ii) is satisfied by taking the minimum over An^ \ A^2^.. The volumes are connected and we 


(46) 


'^min 




1 - Af 


1 - Af ^ 
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This is of the form given in (TT^ . so we choose 


(47) 


Si, ■= 


Cr, 


C„ 




\ 1-^: 




for the smallest value of £j > 2 such that ej. < For each j 1 the martingale method in the Xj direction 
gives a bound of the form 

(1 


(48) {bj+i,- ■■,bd)\ Al^J{bj+u ■■■,bd)) 

O-j 


An, 


Martingale Method in the xi Direction: For the last application of the martingale method to 
^'b 2 +e 2 (^ 3 : ■■■,bd)\ Af'^ (& 3 ,..., 6 d), we use the sequence 

d 

A^^\b 2 ,...,bd) = ... ,bd) \ A^^J {bs,... ,bd)^ n {v x < n + '^{^ - Vj){-bj)} 


t=2 


= {xGZ‘^ : 0 < 


mi 


< n + ^(— - Vj){-bj), for j = 2,...,d, bj < xj < bj+£j}. 


i=2 


mi 


The form of the equations from this set are the same as those from A' where we have replaced Li with n, and 
L with —bj. As such, the upper and lower bounds from C'(A') hold for C{An\b 2 , • ■ ■, bd)) after making the 
appropriate substitutions. We will use this fact when discussing Condition (iii). We note that this volume 
has length n in the the xi direction when Xj = bj for all j 1, and as such the bounds on xi are well 
defined. 

For Condition (ii) we need to consider the minimum of the spectral gaps of two finite families of volumes. 
For n > £i, the volumes have the form 

AW( 62 ,...,M\Ai' 2 ,^( 62 ,...,M 

d 

:= {x ■■ for j = 2,..., d, bj < Xj < bj + £j,n — £i < v ■ x + ^(—- — Vj)bj < n}. 

j=2 

These volumes are isomorphic for any choice of n, 62 , ■ • ■, bd since the slanted boundary generated by v only 
has nonzero components 1 and —1. We denote by P{£) a representative of these volumes. For n = £ 1 , 
A-nlii = = 0 by convention. In this case we are left to consider 

a«(& 2 ,...,M 

d 

< ^1 + - ^>j)i-bj), for j = 2 ,..., d, b, < Xj < bj + £j}. 

j=2 

We denote the resulting finite family of volumes as T(^. Condition (ii) is satisfied by taking a minimum of 
spectral gaps over the union of both sets 

(49) 


:= {x G 0 < vx 

nil 


:= min 7 (A). 
AG{P( 7 )}ur( 7 ) 


For Condition (iii), the volumes are connected and we may apply Lemma 
recall that 

C(aW AAi^)) ^ ^ 
C(AWi\AiV^A“ ■ 

Using the bounds on C(a 1^^) as previously stated, we bound the norm by 

.(II. .(II _ 

(50) 


To simplify the bound. 


C(An ) ^min 
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which decays exponentially as Ai >1. We choose 
(51) := max(Ai, 

for the smallest value of > 2 such that 

Combining the Bounds: As in Case 3a, we iterate the bounds in each direction. The bound on 7 (Al) 
is 


(52) 


d 

7 (Al) > min 7 (A) TT 
AG{p(^)}ur(^) fj. 




> 0 . 


Using Theorem 13.11 this is a positive lower bound for j{D). This completes the proof of the lower bounds 
for all cases. 

Sketch of Proof for Corollary 12.31 Suppose log A ^ 0. Choose a unit vector m such that log A ^ 
±11 log A||m and for which there exists j € {1,..., d} such that Aj ^ 1 and ^ 0 Then, the conditions of 
Case 3a are satisfied. Since the PVBS models are translation invariant, the sequence of volumes A^, from 
Case 3a may be translated in such a way that A^ ^ Therefore, by Theorem 13.II 

7(Z'^) > liminf7(AL) > 0. 


6 . The Upper Bound on the Spectral Gap 


Recall that is one-dimensional and spanned by It follows that any one particle state in the GNS 
Hilbert space is orthogonal to the ground state space. By the variational principle the energy of any 
state in'H.D that is orthogonal to Qd is therefore an upper bound for the spectral gap. That is, if 'I' ± C/£), 
then 


l{D)< 


II^P 


We prove Theorem 12.11 by applying the variational principle to the sequence of one particle states, of 
the form 


where ttq is the GNS representation and A^ is defined by 

Al := {x € ■. 0 < m ■ X < L, for j = 2,..., d, —L < Xj < L}. 


Note that 7 ro(cr^)Ho is the product state describing a particle localized at the site x, and that Al Z' D. For 
any observable A with finite support, the GNS representative 7 ro(A) acts on the product vector in the 
canonical way. From now on, we will drop ttq in the notation. 

Theorem 12.ll is proved by establishing that for all d > 2, A £ (0, 00 )'^, and unit vectors m £ R'^, 

(53) limsup ^^ f ^ ||logA|||sin( 6 >)|, 

ll^^lr c(m)c(A)2 

where 0 is the angle between log A and the outward normal —m of the hyperplane and c(*) is defined as 
in 0 . We restrict to the case that 9 £ (— 7 r/ 2 , 7 r/ 2 ) as for all other values of 9, the lower bounds derived in 
Sections m and [S] show that the gap does not vanish, and we are interested in the behavior near 9 = 0. 


Proof. The expression for the upper bound given in the theorem is invariant under the permutations and 
reflections. We use the coordinate transformations discussed in Remark 12.41 to assume rrij > 0 for all j. Since 
we assume m • log A < 0, there must be at least one value of the index j such that nij > 0 and Xj £ (0,1). 
Without loss of generality we assume this holds for j = 1. 

We prove a lower bound on and an upoerbound on to obtain (1551) . Let A^ = 

for j ^ 1. For a fixed X 2 ,... ,Xd in A^, we denote the minimum value of xi by a{x) = — Yl'j =2 
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Figure 4. A typical in d = 2. An edge between a site on the boundary of A^ and a 
site outside of A^, indicates a nearest neighbor interaction that contributes nonzero energy 


to 


where r < 1 is a remainder term that depends on X 2 , ■ ■ ■ ,Xd- Since Ai < 1, it follows that > Af, and we 
find 


(54) 


L L a{x)-\-L — l 

= E E E 

Xd——L X2 — — L xi—a{x) 

E E 

Xd — — L X2 — — L aii=0 



We note that C{Kl) = diverges as Al Z' D since each sum J2xj=-L diverges. Appealing to 

Proposition 2.2], shows the lower bound on ||4'^|p explicitly proves that Qd = span{0^}. 

For the upperbound on since 


d 

H = E E ^x,x+ej 

j—1 x,x-\-ej^D 


and is the one particle ground state on Aj;, and the zero particle ground state on D\Al, the only nonzero 
contributions to the energy are given by interactions with support {a;, y} such that x G A^ and y G D \ A^, 
see Figure ID Since the PVBS model is nearest neighbor, it follows that there exists a j such that either 
y = x + Cj or y = X — Cj . Therefore, 

= ( 4 -^, ^ 

\x-y\ = l 

Note that for a site x G Al to have a nearest neighbor outside Al requires that x is on the boundary of 
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Ai. Computing hx,y'^^) for x on the boundary and y = x + Cj or y = x — ej outside of produces 




x^ + 1 


A| + l 


< A 


< A^ 


2x 


We consider each boundary that is connected to D\Al separately. For j > 2 and x G Al such that Xj = L, 
the only nearest neighbor y G D\ A^ is y = x + ej. Similarly, for x G A^ such that Xj = —L, the only 
nearest neighbor y G D \ A^ is y = x — ej. For x G A^ such that xi = a{x) + L — 1, there is a maximum 
of d nearest neighbors y G D \ A^, namely yi = x + ei,... ,yd = x + Cd, see FigurelH Since the interaction 
terms of the PVBS Hamiltonian are non-negative, it follows that 




j=2 xGAl 
Xt—L 


j=2 xGAl 
Xi——L 




xGAl i=i 
xi—a{x)-\-L — l 


(55) 


^ E E 


^E 


x£Al 

Xi =±L 


x£Al 

i:i—a{x)-\-L — l 


For a fixed value of j = 2,... d, we have 


E E E 

x^Al Xj—:kL Xd——L 

Xi=^L 


L a{x)+L — l 

■ ^ Ar= E 

X 2 — — L xi—a{x) 


(56) 

where we use Xf' < 1. Similarly, 


<(Af+ A-2^) 


n E 

yk^l,j Xk = -L 


l-Xf^ 

1-Xl 


L L 

y2xd _ ^ ^2ai2^2(a(a:)+L-l) 

Xd——L X2 — — L 


x£Al 

xi—a{x)-\-L — l 


(57) 


S n E >i“ 


V/c—2 X——L 


We replace (l56)l and (|57)l into (155]) to obtain the final upper bound for Using the lower bound 

from (I54p . we find that the overall energy contribution is bounded above by 


(58) 




^A2(i-1)(J_A2) 


1-Af 


X2xj\ 

i=2 yl^xj^-L'^j ) 


The first term is exponentially small in L, and tends to zero. For each j > 1, if Aj = 1, then the 
corresponding term in the sum equals 2j{2L + 1) and converges to zero as L goes to inhnty. If Xj ^ 1, 
treating the cases Aj < 1 and A^ > 1 separately at taking the limit L —^ oo we find 


E Z/ T 2X7 

x,=-L \ 


1 - min(Aj, ) . 
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rrij 

Using the definition Xj = XjX^ ^ can rewrite min(Aj, X~^)‘^ = e ^ mi Note that in the 

case that Xj = 1, we have 0 = 1 — miii(Aj, Therefore, by the variational principle 

^ rrii 

7(^) < Ar" ^ 1 - ^' 

2 ^ 

— \2 E I — rrij log Ai + mi log Xj \ 

j=2 

where the last inequality holds using the bound 1 — e““ < x. In order to make this estimate independent of 
the choice of Xi, we define c{v) = min{|?;j | : Vj ^ 0} for u G and bound mi and Af below by c{m) and 
c(A), respectively. 

The expression —m^TogAi TmilogAj can be interpreted as (logA,/j), where fj = —mjCi + miCj, 
j = 2,...,d. Let P denote the orthogonal projection onto the bounding hyperplane perpendicular to m. 
Clearly {m, fj) = 0, hence Pfj = fj. We have ||PlogA|| = || log A||| sin0|. Also, note \\fj\\ < ||m|| = 1. With 
this, we can make the following estimate: 

d d d d 

= E = E ^EII^^°sA|| < (d-1)11 log A|||sin6»|, 

j=2 j=2 j=2 j=2 

where 9 denotes the angle between —m and log A. 

7(g^)< J |logA|||sin(d)|. 

c(m)c(A)^ 

□ 


We now prove of Lemma 13.31 


Appendix 


Proof. We bound the norm of the product of projections ||G^"+i\^"+i-'£’„||. It sufficient to consider 
GA'*+i\A"+i-^ acting on the range of En, that is C/a„ fl A vector in this subspace has the form 

(59) T = ® ^ ® CTi4'An+i\A„ ^ ^ ^i^a„+aa„ 

irGAn+i\A„ a;GAn,+i\A„, 


where bo, Ox, G C, Tq , diA a,re the zero and one particle ground states in finite set A, and (t,),'Lq 
is the product state with a single particle at site x G A„+i \ A„. 

By construction, the vector is orthogonal to . For to be perpendicular to the single particle 

ground state d/A^+i |;]^g coefficients must satisfy the following orthogonality condition: 


(60) 


bo = 


-I 


E 

xGAn4-i\ATi 


UxAP 


Since A„+i\ A„+i_£ is connected, the ground state space is two dimensional. The projection G'An+iVA^+i.j, 
is onto the zero particle ground state vj/A^+^^An+i-^ single particle ground state, v[/An+i\A„+i_^ ^ 

drop the notation A„+i \ Kn+i-e. from the ground state vectors when the volume is understood. The 
expression GA^+i^An+i-^il/ has four nonzero terms: 

1. The projection onto a vector with no particles in A„+i_^. 

(61) |d'o)('ko|d' = 0 ,pWaAa„a.-. 
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2. The single particle projection onto the bo term. The nonzero terms have a particle in A„ \ Kn+i-f. 


(62) 




A„+i\A„ 


= 6o 


3. The single particle projection onto the 


(63) |vI/i)(M/i| ^ a.|l).= 


aieAn+i\A„ 




_ C'(A„ \ Kn+i_i) _^ ^A„+i\A„+i_f 

v'C(A„)C(A„+i \ 

terms. The nonzero terms have a particle in A„_|_i \ A„: 

n+l \ ^n+l-t) a:GA„+i\A„ 


4. The single particle projection onto the terms. The nonzero terms have a particle in A„+i \ A„ and 
another in A„+i_^: 


(64) ^ 6,^'^0|1), 

a:eAn,+i\An 


\/ C{A.n+l-l) 

v'C(A„)C(A„+i\A„+i_,) 


E 

rr’C. \ . « \ A 




a+l-l, 


^^A„+i\A„+i_, 


The inner product of the the sum of these vectors gives the (square) of the norm; 


||(^A„+i\A„+i_,^||2 


|, ,2 ^{^n+l-t) _ C{Kn+l-t) _ 

' C(A„) ^C(A„)C(A„+i\A„+i_,) 


6oC(A„ \ A„_|_i_^) ^ 1 ^ 

■\/C'(A„)C'(A„+i \ A„_|_i_f) yj C(A„+i \ A„+i_^) 


^ b.X 

A„+i\A„ 


Applying the orthogonality condition for |5op and combining the Qx terms simplifies 


(^(An+i-f) 

C(A„)2 


y] a,A" 

2 

+ 

A„+i\A„ 



C(A„)-C(A„\A„+i,) 


^ ax^ 
A„+i\A„ 

a^A^ 

A„+i\A„ 

y] ax\^ 

A„+i\A„ 

y] a,A^ 

A„+i\A„ 

y^ “s^A^ 

A„+i\A„ 


C(A„)\/^^(A„+i\ -^n+1—4 -i\A 


C^h-n+l-tf 


y] axA^ 


C{h-n+l-i) 

L C(A„)2 + C(A„)2C(A„+1 \ -^n+1—£) J 


^(A^+i.,) (C(A ri+l \ A„+i_f ) + C(A„+i_,)) 


C(A„)2C(A„+i \ K+i-t) 


C(A„+i,)C(A„+i) 


LC(A„)2C(A„+i \ An+1—£) J 


C{K+l-t) (C(An+l \ An) + C(A„)) 
C(A„)2C(A„+i \ ■^n+1 —£) 

C(A„+i^) 


C(^h.n)C{^Kn+l \ X.ji+1—i) 


C[K+i\K) 

C(A„) 
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We combine terms and apply Cauchy-Schwarz for an upper bound: 


IIGvPf = 




< 


C(A„)C(A„+i \ 

G(A„+i_i)C'(A„+i \ A„) 
C(A„)C(A„+i \ 

-^n+1—f) 


A„ + AAn 


^ a.y 

A„+1\A„ 


+ 


C{kn) 


A„+i\A„ A„+i\A„ 


IE 


A„+i\A„ 


G(A„) 


_ G(A„-|-1_^)C'(A,^-|-1 \ Kn) ||, t ,||2 
“ C(A„)C(A„+1 \ -^n+1—f) 

We have equality when we choose Ux = bx = . Therefore, 


II ^ p |j2 ^(-^n+1—^)^(-^n+l \ -^n) 

"" " C(A,)C(An+i\An+i-£)’ 
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